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We discuss hadron production in heavy ion collisions at RHIC. We argue that hadrons at transverse 
momenta Pt < 5 GeV are formed by recombination of partons from the dense parton phase created 
in central collisions at RHIC. We provide a theoretical description of the recombination process 
for Pt > 2 GeV. Below Pt = 2 GeV our results smoothly match a purely statistical description. 
At high transverse momentum hadron production is well described in the language of perturbative 
QCD by the fragmentation of partons. We give numerical results for a variety of hadron spectra, 
ratios and nuclear suppression factors. We also discuss the anisotropic flow V2 and give results based 
on a flow in the parton phase. Our results are consistent with the existence of a parton phase at 
RHIC hadronizing at a temperature of 175 MeV and a radial flow velocity of 0.55c. 



I. INTRODUCTION 

Recent data from the Relativistic Heavy Ion Collider 
(RHIC) have shown a strong nuclear suppression of the 
pion yield at transverse momenta larger than 2 GeV /c in 
central Au + Au collisions, compared to p+p interactions 
lj. This is widely seen as the experimental confirma- 
tion of jet quenching, the phenomenon that high energy 
partons lose energy when they travel through the hot 
medium created in a heavy ion collision 0, ff Q > entail- 
ing a suppression of intermediate and high Pt hadrons. 

However, the experiments at RHIC have provided new 
puzzles. The amount of suppression seems to depend on 
the hadron species. In fact, in the production of pro- 
tons and antiprotons between 2 and 4 GeV/c the sup- 
pression seems to be completely absent. Generally, pions 
and kaons appear to suffer from a strong energy loss while 
baryons and antibaryons do not. Two stunning experi- 
mental facts exemplify this [MISS 111- First, the ratio of 
protons over positively charged pions is equal or above 
one for Pt > 1.5 GeV/c and is approximately constant 
up to 4 GeV/c. Second, the nuclear suppression factor 
Raa below 4 GeV/c is close to one for protons and lamb- 
das, while it is about 0.3 for pions. 

There have been recent attempts to describe the differ- 
ent behavior of baryons and mesons through the existence 
of gluon junctions 9j or alternatively through recom- 
bination as the dominant mechanism of hadronization 
[ToL llT|| . The recombination picture has attracted addi- 
tional attention due to the observation that the elliptic 
flow pattern of different hadron species can be explained 
by a simple recombination mechanism [T3. Il3l 111 . IT^ | . 
The anisotropics v% for the different hadrons are com- 
patible with a universal value of i>2 in the parton phase, 
related to the hadronic flow by factors of two and three 
depending on the number of valence quarks p"fj| . 

In this work we elaborate on the arguments presented 



in [lOj. We will present the formalism as well as new 
numerical results. The competition between recombi- 
nation and fragmentation delays the onset of the per- 
turbative/fragmentation regime to relatively high trans- 
verse momentum of 4-6 GeV/c, depending on the hadron 
species. This is the explanation for several key observa- 
tions at RHIC: 

• the two component form of hadron transverse mo- 
mentum spectra, including an exponential part and 
a power law tail with a transition between 4 and 6 
GeV/c. 

• the very different behavior of the nuclear suppres- 
sion factors Raa of mesons and baryons. 

• the particle dependence of the elliptic flow. 

• the unusually large baryon/meson ratios. 

In addition, our calculation clarifies the range of applica- 
bility of perturbative calculations including energy loss. 

The paper is organized as follows. In the next section 
we explain why fragmentation might not be the domi- 
nant mechanism of hadronization at intermediate trans- 
verse momenta of a few GeV/c in heavy ion collisions. 
We discuss the fundamental principles of recombination 
and fragmentation. In Sec. Ill we present the theoreti- 
cal framework of recombination in more detail and also 
discuss its shortcomings. In particular, we address the 
question of applicability at low transverse momentum. 
In Sec. IV we introduce our parametrization of the par- 
ton spectrum and discuss further calculational details, 
and in Sec. V we present numerical results on spectra, 
hadron ratios, nuclear suppression and elliptic flow. Sec. 
VI summarizes our work. 
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II. FRAGMENTATION VS. RECOMBINATION 

A. Fragmentation of partons 

Inclusive hadron production at sufficiently large mo- 
mentum transfer can be described by perturbative quan- 
tum chromodynamics (pQCD). The invariant cross sec- 
tion for a hadron h with momentum P can be given in 
factorized form 

a Q 

The sum runs over all parton species a and o~ a is the cross 
section for the production of parton a with momentum 
P a = P/z. Thus the parton production cross section 
has to be convoluted with the probability that parton a 
fragments into hadron h. The probabilities D a ^h(z) are 
called fragmentation functions • Like parton distribu- 
tions they are non-perturbative quantities. However they 
are universal and once measured, e.g. in e + e~ annihila- 
tions, they can be used to describe hadron production in 
other hard QCD processes. 

Using Q we can estimate the ratio of protons and pi- 
ons. Taking, e.g., the common parametrization of Kniehl, 
Kramer and Potter (KKP) [lj|, the ratio D a ^ p / Da^^a 
is always smaller than 0.2 for each parton a. This re- 
flects the well known experimental fact that pions are 
much more abundant than protons in the domain where 
pQCD is applicable. The excess of pions over protons 
even holds down to very low Pt, smaller than 1 GeV, 
where perturbative calculations are no longer reliable. In 
that domain one can argue that the difference in mass, 
M p 3> M w , lays a huge penalty on proton production. 
The small value of the p/ir° ratio predicted by these cal- 
culations over the entire range of Pt is the reason why 
the ratio p/tt° ~ 1 measured at RHIC is so surprising. 

It has been suggested that the fragmentation functions 
D a ^h( z ) can be altered by the environment poll2ll|. The 
energy loss of the propagating parton in the surrounding 
medium leads, in first approximation, to a rescaling of the 
variable z. This would affect all produced hadrons in the 
same way, and thus cannot explain the observations at 
RHIC. In a picture with perturbative hadron production 
and jet quenching alone, the different behavior of hadrons 
can not be described by one consistent sets of energy 
loss parameters. To save the validity of the purely per- 
turbative approach species dependent non-perturbative 
contributions to the fragmentation functions have to be 
introduced ad hoc to explain the data |22j . 

Perturbative hadron production consists of three steps: 
production of a parton in a hard scattering, propagation 
and interaction with a medium, and finally hadronization 
of the parton. Only modifications in our understanding 
of hadronization are able to provide an explanation of the 
experimental observations, since the other steps are blind 
to the hadron species that will eventually be created. 



B. From fragmentation to recombination 

For the production of a hadron with momentum P via 
fragmentation we need to start with a parton with mo- 
mentum P/z > P. The fragmentation functions favor 
very small values of z, i.e. the situation where the en- 
ergy of the fragmenting parton is not concentrated in 
one hadron. On the other hand, the transverse momen- 
tum spectrum of partons is steeply falling with Pt- This 
makes it clear that fragmentation is a rather inefficient 
mechanism for the production of high Pt hadrons, since 
it has to overcome the limited availability of partons at 
even higher transverse momentum. As a result, the av- 
erage (z) is larger than what is expected from the shape 
of the fragmentation functions. For pion production, (z) 
is about 0.6 for the production from a valence quark, 
0.4 for a sea quark and 0.5 for a gluon in the range 
2 GeV/c < P T < 10 GeV/c for leading order KKP frag- 
mentation functions. 

An outgoing high energy parton is not a color sin- 
glet and will therefore have a color string attached. The 
breaking of the string will initiate the creation of quark 
antiquark pairs until there is a entire jet of partons, which 
have to share the energy of the initial parton. They will 
finally turn into many hadrons. The creation of several 
hadrons from one fragmenting parton is the reason why 
fragmentation functions prefer small values of z. If phase 
space is already filled with partons, a single parton de- 
scription might not be valid anymore. Instead one would 
have to introduce higher twist (multiple parton) fragmen- 
tation functions. In the most extreme case, if partons are 
abundant in phase space, they might simply recombine 
into hadrons. This means that a u and a d quark that 
are "close" to each other in phase space can bind to- 
gether to form a ir + . The scale of being close will be 
set by the width of the pion wave function. In this sce- 
nario the total pion momentum will be just the sum of 
the individual quark momenta. We immediately notice 
that this recombination mechanism is very efficient for 
steeply falling spectra: in order to produce a 5 GeV pion 
we can start with two quarks having (on average) about 
2.5 GeV/c transverse momentum and each being there- 
fore far more abundant (on average) than a 10 GeV/c 
parton that could produce the pion via fragmentation. 
Of course the recombining partons must be close in phase 
space, i.e. recombination will be suppressed if the phase 
space density is low. 

Recombination can be interpreted as the most "exclu- 
sive" form of hadronization, the endpoint of a hypotheti- 
cal resummation of fragmentation processes to arbitrary 
twist. We cannot achieve a quantitative understanding 
of this at the moment. However, we can try to find an 
effective description which can be tested against observ- 
able consequences. In this work we will advocate a sim- 
ple model for recombination and compare it with single 
parton fragmentation. These two mechanisms of hadron 
production compete differently depending on the phase 
space density of partons. From the above we understand 



3 



that the competition between fragmentation and recom- 
bination is dominated by the slope and the absolute value 
of the phase space distribution of partons. Below we show 
that recombination always wins over fragmentation for an 
exponentially falling parton spectrum, but that fragmen- 
tation takes over if the spectrum has the form of a power 
law, as it is provided by pQCD. We will apply this insight 
to hadron production in relativistic heavy ion collisions 
at midrapidity and transverse momenta of a few GeV/c 
where we expect a densely populated phase space. For 
the recombination of three quarks into a proton the mo- 
menta of three partons have to be added up, but only 
two momenta in the case of a pion. Assuming an expo- 
nential parton spectrum this implies for a proton a distri- 
bution ~ [exp(-Pr/3)] 3 and for pions — [exp(-PT/2)] 2 , 
predicting a constant p/n + ratio where the value is de- 
termined by simple counting of quantum numbers |10| . 
We will show that some of the surprising experimental 
results from RHIC can be explained in this way. 

C. The recombination concept 

The idea of quark recombination was proposed long 
ago to describe hadron production in the forward region 
of p + p collisions 23]. This was later justified by the 
discovery of the leading particle effect, the phenomenon 
that, in the forward direction of a beam of hadrons col- 
liding with a target, the production of hadrons sharing 
valence quarks with the beam hadrons are favored. E.g. 
in the Fermilab E791 hxed target experiment with a 500 
GeV 7r~ beam [24j the asymmetry 

da D -/dx F - da D +/dx F 

a ( XF > = i u — z~a u — ( 2 ) 

dajj- jdxp + d(T£>+ /dXF 

between D~ and D + mesons grows nearly to unity when 
the Feynman variable xf, measuring the longitudinal 
momentum relative to the beam momentum, approaches 
1. Fragmentation would predict this asymmetry to be 
very close to zero. However, recombination of the c quark 
from a cc pair produced in a hard interaction with a d 
valence quark from the tt~ , propagating in forward direc- 
tion with large momentum, is highly favored compared 
to the recombination of the c with a d which is only a sea 
quark of the 7r~. This leads to the enhancement of D~ 
over D + mesons in the forward region. 

The leading particle effect is a clear signature for the 
existence of recombination as a hadronization mechanism 
and has been addressed in several publications recently 
[25L |2(j . In this case recombination is favored over frag- 
mentation only in a certain kinematic situation (the very 
forward direction), which is a only a small fraction of 
phase space. 

In central heavy ion collisions many more partons are 
produced than in collisions of single hadrons. The idea 
that recombination may then be important for a wide 
range of rapidities - and at least up to moderate trans- 
verse momenta - was advocated before |2j], |28|, |29J ■ How- 



ever, it was only recently that RHIC data indicated that 
recombination could indeed be a valid approach up to 
surprisingly high transverse momenta of a few GeV jc. 

Charm and heavy hadron production have the advan- 
tage that the heavy quark mass provides a large scale 
that permits a more rigorous treatment of the recom- 
bination process |26| . The description of recombination 
into pions and protons seems to be theoretically less rig- 
orous. However, a simple counting of quantum numbers 
in a picture where the structure of hadrons is dominated 
by their valence quarks often provides surprisingly good 
results. This has been pointed out for particle spectra 
and ratios [H |H and for elliptic flow [H [H E3 . 
Most of the work so far has stayed on a qualitative level 
without quantitative predictions. Recombination of D 
mesons in heavy ion collisions has been investigated in 

M 

We will argue below that the counting of quantum 
numbers is a good description of the recombination pro- 
cess for intermediate momenta. The formalism will set 
the stage to obtain quantitative results in this regime by 
recombining quarks from a possibly thermalized phase. 
We know that the parton phase will not behave like a 
perturbative plasma near the hadronization point |33| . 
Instead quarks at hadronization will be effective degrees 
of freedom exhibiting a mass and gluons will disappear 
as dynamical degrees of freedom. We will assume here 
that the effective quarks behave like constituent quarks 
and that there are no dynamical gluons. 

This is different from the work of Greco et al. [Ill ] who 
suggested to recombine one perturbative quark with ther- 
mal quarks, leading to an additional contribution at the 
transition region between the pure thermal phase domi- 
nating below 5 GeV/c and the pure fragmentation regime 
dominating above 5 GeV jc. A similar form of this pick- 
up reaction of a perturbative parton was recently pro- 
posed in the context of the sphaleron model 34] . 

One can also attempt to extend the recombination con- 
cept to low Pt, however, the theoretical situation is much 
more ambiguous there. The main reasons are that the 
simple counting of quantum numbers violates energy and 
entropy conservation at low Pt, where the bulk of the 
hadron yield resides. Nevertheless, once recombination 
is recognized to be the dominant hadronization mecha- 
nism at intermediate Pt from 2 to 5 GeV/c, it is quite 
reasonable to expect that this mechanism extends down 
to very small transverse momentum, where quarks are 
even more abundant. However the theoretical descrip- 
tion will only be on solid ground once the problems of 
energy and entropy conservation are addressed properly. 



D. A non-relativistic model 

For a first estimate, let us consider a simple static 
model for the recombination process. We start with a sys- 
tem of quarks and antiquarks which are homogeneously 
distributed in a hxed three-dimensional volume V with 
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phase space distributions w(p), so that the number of 
quarks or antiquarks with a particular set a of quantum 
numbers (color, spin, isospin) is 



N a = V 



d 3 P 



Wa(p) 



(3) 



Hadronization is assumed to occur instantaneously 
throughout the volume. The distributions w(p) are sup- 
posed not to change during the hadronization process, 
i.e. the quarks are assumed to be quasi- free. 

The spatial wave functions for a two particle 
quark/antiquark state with momenta pi and P2 and a 
meson bound state with momentum P are 



(x\q, Pl p 2 ) = y-V(pi*i+P2x 2 ) 
{x\M,-p)=V- 1 / 2 e ip - R i PM (y) 



(4) 
(5) 



respectively, with the center of mass and relative coor- 
dinates R = (xi + X2)/2 and y = xi — X2 for the two 
quarks in the meson system. To keep our notation simple 
we omit the proper antisymmetrization of multi fermion 
states, since all combinatorial factors will cancel in the 
final result. The internal meson wave functions is nor- 
malized as 



d 3 y\<p M (y)f 



1. 



(6) 



The overlap amplitude is given by 
(q,pip 2 |tt,P) 



y3/2 



■<5 3 (P-pi-p 2 )£ M (q). (7) 



Here, we have introduced the relative momentum q = 
(Pi — P2)/2, conjugate to y, and <p M (q) is the Fourier 
transformed wave function. The squared amplitude is 



|<g,piPi|M,P) 



V 2 



<J 3 (P-pi-p 2 )|£M(q)| • (8) 



We conclude that the total number of pions found in 
the quark/antiquark distribution is 



N M = C M V £ 



d 3 P d 3 Pl d 3 p 2 
(2tt) 3 (2tt) 3 (2tt) 3 



x W (p 1 ) W (p 2 )|(g,p 1 p 1 |M,P)| (9) 

with a degeneracy factor Cm- The momentum distribu- 
tion of the pions is given by 



dN 



A I 



d 3 P 



V 



d 3 q 

(2tt) 3 J (2tt) 3 



x w 



q w 



-q ) hMq) ( 10 ) 



From above equation the spectra of mesons can be cal- 
culated for given quark distributions. This will require 
knowledge of ipM- The wave function has some width 
Am, and we assume that it drops rapidly for |q| > Am- 
Let us study the kinematic region where P ^> Am- The 
integral over the relative momentum q is dominated by 
values |q| ~ Am and thus we can assume that |q| <C |P|. 
We apply a Taylor expansion 



/P \/P \ / P \ 2 " / P \ /P\ /P\ /P\ 

{2 +q ) W {2- q ) =W {2) W {j) d ^ W {2)~ dlW {j) d H2) 



0(q 3 



(11) 



r 



where the first order term in the expansion vanishes. 

From the lowest order term in the expansion we get a 
contribution to the meson spectrum which is independent 
of the shape of the wave function. Only the normalization 
© enters and leads to a universal term 



(12) 



The second order term (like all higher order terms) in the 
expansion generates a correction depending on the shape 
of the wave function. For the sake of simplicity, let us 
assume a Gaussian shape 



<£iw(q) =A4re' 



-q 2 /2A? 



The normalization factor is determined by 10 as N M — 
(2y/ir I Am) 3 . For this example we have the second order 
correction 



V 



A 2 



(2tt) 3 



w(P/2)Aw(P/2) - (Vw(P/2)) 2 . (14) 



We want to emphasize once more that in the limit 
P Am the exact shape of the wave function is not 
important. In a relativistic framework this statement is 
softened by the fact that we consider the hadron forma- 
tion in a boosted frame, where P is large. Therefore Am, 
which is of order Aqcd in the rest frame of the hadron, 
will be dilated. 

As an example, let us assume an exponential parton 



(13) distribution of the form w(p) = e 



-P/T. 



The meson spec- 
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trum at large P is then given by 



d?P 



Ctt 



V 



2*) 



-P/T 



OA 2 1 



(15) 



using l|12 |l .l[Tl jl . The second order term introduces a 
power correction of order A M /TP to the universal re- 
sult, depending both on the width of the wave function 
and the slope of the parton distribution. 

For nucleons we start with three quarks at coordinates 
Xi . We introduce center of mass and relative coordinates 
R = (xi + x 2 + x 3 )/3, y = (xj. + x 2 )/2 - x 3 and z = 
x a - x 2 . 

The overlap amplitude between a 3 quark state and a 
baryon with momentum P is 



(q,pip 2 p 3 \B,P) 



(27T) 

V 2 



-^ 3 (P-Pi-P2-P3)^s(q, s). 

(16) 



Here <^s(q, s) is the baryon wave function in momentum 
space, depending on the relative momenta q and s con- 
jugate to y and z respectively. 

Hence the baryon distribution is given by 



dN B n V 



d 3 q d 3 s | „ . .,2 
Vs(q,s) 



(2tt) 3 J (2tt) 3 (2tt) 3 



x w 



/P + q 

V 3 2 



(17) 



Cs is the appropriate degeneracy factor. In the region 
where the nucleon momentum P is large compared to 
the intrinsic width of the wave function, we can again 
expand the product of the quark distribution functions. 
The leading term is universal and contributes 



V 



(2tt) : 



-i«(P/3) ; 



(18) 



to the nucleon distribution. The second order term pro- 
vides a correction 



v 1 1 , 

(2^)3 2 1 m 



h 2 



;(P/3) 2 Aw(P/3) - w(P/3)(Viu(P/3)) : 



(19) 



assuming a normalized Gaussian wave function 

<^ S (q,s) =W i3e -<l 2 /2A 2 Ble -s 2 /2A 2 B2 . (20) 



E. First conclusions 

Let us summarize and analyze our first results. The 
transverse spectrum of mesons from recombination is 
proportional to Cmw 2 (Pt/2) whereas fragmentation 



would generate a distribution proportional to D(z) ® 
w(Pt/z). For an exponential parton spectrum w — 
e -P T /T ratio of recombination to fragmentation is 



R _ C M -- 
F (D) e 



(21) 



where (D) and (z) < 1 are average values of the frag- 
mentation function and the scaling variable. Obviously, 
for large Pr one always gets R/F > 1. In other words, 
recombination always wins over fragmentation at "high" 
Pt for an exponential parton spectrum. The same is true 
for baryons as well as mesons. 

Now let us consider a parton distribution given by a 
power law spectrum w = A(Pt / ' n)~ a with a scale \x and 
a > 0. Then the ratio of recombination over fragmenta- 
tion is 



R CmA /_4_ 

F (D) \(z) 



Pt 
Q 



(22) 



and fragmentation ultimately dominates at high Pt- 
Again, this holds both for mesons and baryons. This im- 
plies that fragmentation from an exponential spectrum 
and recombination from a power law spectrum are sup- 
pressed. We will thus omit these contributions in our 
work. 

One might ask whether these considerations are still 
valid in the case of more realistic formulation of recom- 
bination. It turns out that the basic formulae obtained 
above are still valid in a relativistic description. Devia- 
tions are less than 20% for Pt > 2 GeV/c. 

Given an exponential parton spectrum we note that 
recombination yields a a constant baryon-to-meson ra- 
tio. The ratio is then only determined by the degeneracy 
factors 



dN§ 

In* 



Cb_ 

Cm 



(23) 



Just counting the hadron degeneracies, the direct p/7r° 
ratio (neglecting protons and pions from secondary 
hadronic decays) would be ~ 2, in contrast to ~ 0.2 from 
fragmentation in pQCD. We will later see that finite mass 
effects and superposition with the fragmentation process 
will bring down the p/~K° ratio from 2 to approximately 1 
in the range between 2 and 4 GeV/c transverse momen- 
tum. 



III. THE RECOMBINATION FORMALISM 

In this section, we turn to a better description of re- 
combination. This will require a more realistic model of 
the parton phase including longitudinal and transverse 
expansion as well as an improved space-time picture. 

Let us consider a system of quarks and antiquarks 
evolving in Minkowski space. We choose a spacelike hy- 
persurface X on which recombination of these partons 
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into hadrons occurs. In the simplest scenario, that could 
be just a slice of Minkowski space with fixed time to, 
leading back to the case we described in the previous 
subsection. 

It has been discussed in the literature [35ll3^| . how the 
freeze-out can be smeared around the hypersurface to ac- 
count for a finite hadronization time. However, RHIC ex- 
periments suggest a very rapid freeze-out. The measured 
two-particle correlation functions are consistent with an 
extremely short emission time in the local rest frame, 
suggesting a sudden transition after which the individual 
hadrons interact only rarely |3?|. 

This can be understood by the fact that the hadroniza- 
tion time in the laboratory frame is Lorentz contracted. 
The Lorentz factors -fr are 3.35 for a 3 GeV/c proton and 
22 for a pion of the same Pt- The hadronization time 
— even if it is more than 1 fm/c in the rest frame of the 
hadron — will be experienced as considerably shorter by 
fast hadronizing particles. 



A. Wigner function formalism 



nuclear collisions Here, we will provide a deriva- 

tion for the recombination of a quark antiquark pair into 
a meson. The generalization to a three quark system 
recombining into a baryon is straightforward. 

By introducing the density matrix p for the system of 
partons, the number of quark-antiquark states that we 
will interpret as mesons is given by 



/d 3 P 
(M;P|p a6 |M;P) 



(24) 



Here |M;P) is a meson state with momentum P and 
the sum is over all combinations of quantum numbers 
- flavor, helicity and color — of valence partons that 
contribute to the given meson M. We insert complete 
sets of coordinates 



M 



d 3 P 



tfh^ficFhd 6 ^ (M;P| fi,f 2 



(2tt) 3 

x (ri,fa|p|fi,^)<fi,^|M;P). (25) 



Recombination has already been considered before in 
a covariant form utilizing Wigner functions for the pro- 
cess of baryons coalescing into light nuclei and clusters in 



and change the variables to ri )2 = (?i,2 + ?i 2)/^ an d 
r i.2 = £1,2— ri 2- We define the 2-parton Wigner function 
W /! a6 (r 1 ,r 2 ;pi,p 2 ) as 



J 



r i r 2 



d 3 pi d 3 p 2 
(2tt) 3 (2tt) 3 



e- ipiri e- l P 2r ^ afc (ri,r 2 ;pi,p2) 



(26) 



and introduce a notation for the meson wave function 



e -,p.(R + R72)^ M ( r _£. 



(27) 



It is convenient to change coordinates again to 



RW = (rW + rW)/2, 



with conjugated momenta 

P = pi +p 2 , 

q = (pi P2V2. 



(28) 
(29) 



(30) 
(31) 



We arrive at 

N M 



4-7 (27T) 3 (27T) 3 



r r P P 

x V^ fc (R+-,R--;-+q,--q 



xe- (p M \r+- up M i r- - 



(32) 



The integration over R' has been carried out and provides 
the three- momentum conservation P = P. 

We define the Wigner function &Yi °f the meson as 

<(r,q)= J dVe-"'^ A ,m^^). (33) 

Then 



dN 



M 



d 3 P 



(34) 
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To evaluate this expression, we have to model the 
Wigner functions of the parton system and of the me- 
son. We assume that the 2-parton Wigner function can 
be factorized into a product of classical one-particle phase 
space distributions w. Furthermore the color and helicity 
states for each flavor will be degenerate. We can there- 
fore replace the sum over quantum numbers a and b by 
a degeneracy factor Cm- 

We introduce the following simplifications: the spatial 
width Ar of the hadron wave function, translating into 
a width of will be small compared to the nuclear 
size of the system at hadronization. The phase space 
distributions w(r; p) are steeply falling functions of the 
momentum p, but vary much less with the spatial co- 
ordinate r within the typical size Ar of a hadron. We 
shall assume that the spatial variation of the phase space 
distribution is small on this scale, replacing R±r/2 with 
R. 

In our derivation we implicitly chose an equal time 
formalism by introducing the spatial coordinate states 
|r) at a fixed time. However, in our result the integration 
over the final state phase space 

d 3 Pd 3 R = d 3 Pd 3 R P ■ u(R) /E (35) 

is manifestly Lorentz invariant. Here E is the energy 
of the four vector P and u(R) is the future oriented 
unit vector orthogonal to the hypersurface defined by 
the hadronization volume. This form can be easily gen- 
eralized for an arbitrary hadronization hypersurface E 
[Hill. We have 

„dN M _ n f d 3 R P ■ u(R) f d 3 q 
E ^P- CM ] (2ir) 3 J (2^ 

x w a (R; I - q) *X(q) w b (r, | + qV (36) 
I 



Here we have rewritten the spatially integrated Wigner 
function of the meson in terms of light cone coordinates 
for quark a and approximated it by a squared light cone 
wave function of the meson 

<f> M (x,k x ) = |<Mz,k_L)| 2 . (40) 

The light cone wave functions we introduce here do 
not necessarily coincide with the light cone wave func- 
tions < &lc( ;e 7 fx) used for exclusive processes in QCD 
[40| . There, a hadron is decomposed into a series of Fock 
states of perturbative partons, starting from the valence 



a, b now only denote the flavors of the valence quarks in 
meson M and 

*X(q) = J d 3 r<(r,q) (37) 

is the spatially integrated Wigner function of the meson. 



B. Local light cone coordinates 

The structure of hadrons is best known in the infinite 
momentum frame which is described in light cone coordi- 
nates. If we let the hadron momentum P define the z axis 
of the hadron light cone (HLC) frame, we can introduce 
the light cone coordinates e.g. q + , q± for the relative 
momentum q. Note that we denote transverse momenta 
in the HLC frame i.e. the component orthogonal to P 
— by the label _L, but transverse momenta in the center 
of mass (CM) frame of the heavy ion collision — orthog- 
onal to the beam axis — by the label T. The labels + 
and — always refer to light cone coordinates in the HLC 
frame for a given P. We fix the HLC frame by a simple 
rotation from the CM frame, i.e. P + = (E + |P|)/\/2. 
We reintroduce the momentum k — P/2 — q of parton a 
in the meson. Assuming a mass shell condition k 2 = m 2 
with arbitrary but fixed virtuality m a , we can rewrite the 
integral 

1.0 

d 3 k = d 4 k 2k°5(k 2 - m 2 a ) = dk + d 2 k ± — (38) 

in HLC coordinates. m a will be of order Aqcd, or more 
precisely, of the order of a constituent quark mass. In the 
HLC frame, where we assume that formally P + — ► oo, we 
parametrize k + = xP + with < x < 1. Since k~ <C k + , 
we have k°/k + w 1/V2. We end up with 



(39) 

I 

structure. In terms of this expansion we know that the 
valence Fock state has only a small contribution at scales 
fj, > 1 GeV. fi is given by the momentum transfer in the 
hard reaction which is described by perturbative QCD 
(pQCD). However, because the momentum transfer in a 
hard exclusive reaction has to be spread over all partons 
in the hadron state, the exclusive process itself acts like 
a filter, weighting the lower Fock states more strongly. 
In other words, the contributions of higher Fock states, 
though more likely in the wave function, are generally 
suppressed by inverse powers of the momentum trans- 



„dN M n fd 3 RP-u(R) fdxP+d 2 ^ . + . , + , 

1FF = M J (2tt)3 J V2(2tt)3 M R -> xP+ >*±) $M(x,-k±)w b (R;(l-x)P+,--k x ), 



8 



fer (higher twist). This usually permits a fairly good 
description of hard exclusive processes in terms of the 
lowest Fock state. 

As we will discuss below, we expect the parton spec- 
trum of a heavy ion collision at freeze out to be composed 
of an exponential part at small transverse momentum 
and a power law tail, given by pQCD, at high transverse 
momentum. We will study recombination in the pQCD 
domain in a forthcoming publication. 

Here, we want to focus on recombination from the 
exponential part of the parton spectrum. It will be 
given by a slope 1/T* with a temperature-like param- 
eter T*. T* also sets the scale for the typical momen- 
tum transfer in the parton medium before hadroniza- 
tion. If T* is an effe ctive blue shifted temperature 
T* = \J {1 + /9)/(l — (3)T in an expanding medium with 
physical temperature T and flow velocity f3, then the typ- 
ical scale T < T* will even be smaller. Here we use 
the phase transition temperature at zero baryon density 
Tw 175MeV EJ. 

We cannot expect that perturbative QCD will work 
as a description of partons at the phase transition |33j |. 
What are the quanta that recombine? We know that in 
pQCD for decreasing scales the non-perturbative (long- 
range) matrix elements describing hadrons get more "va- 
lence like" , though we cannot seriously extend this study 
to scales below 1 GeV. We will assume here that we 
recombine effective constituent partons, taking into ac- 
count only the valence structure of the hadron. Gluons 
are no dynamic degrees of freedom in this picture, and 
the quarks and antiquarks will have an effective mass. 

This picture is supported by the recent discovery of 
"magical factors" of 2 and 3 in measurements of spectra 
and the elliptic flow of mesons and baryons respectively 
at RHIC [16|. Later we will apply our assumptions to 
partons in the exponential spectrum having as much as 2 
GeV/c of transverse momentum. This might raise doubts 



about the validity of an effective description. However, 
we have to keep in mind that the momentum is in prin- 
ciple meaningless, only the momentum transfer experi- 
enced by a particle in a reaction sets the scale at which 
we resolve its structure. 

From the normalization condition 

(M ; P|M; P') = (2tt) 3 <5(P - P') (41) 

we infer J d 3 ripM(r) = 1 and finally 

/w'*"^' 1 '- 1 - (42) 

This is different from the light cone wave functions in 
exclusive processes which involve a dimensionful quantity 
connected to the weight of the Fock state, e.g. the pion 
decay constant f n in case of the pion. 

We further utilize the fact that we work in a kine- 
matic regime where P + is large compared to all non- 
perturbative quantities. Of course, our main concern 
here is that we do not really know the shape of the wave 
function <f>M- We can choose a factorized ansatz 

$ M (x,k±.) = cf>M{x)Sl(kx) (43) 

with a longitudinal distribution amplitude 4>m(x) and 
a transverse part. We know that the transverse shape 
should be quite narrow, e.g. given by a Gaussian with 
a width Aj_ <C P + ■ Considering hadrons at mid rapid- 
ity, kj_ will mainly be pointing in the longitudinal and 
azimuthal directions in the CM frame, where the varia- 
tion of the parton distributions w is small. This implies 
w(R; xP + , kjj w w(R; xP + ) for typical transverse mo- 
menta kx ~ Aqcd -C P + ■ We can then integrate the 
dependence of the wave function. This leaves us with 



E 



dN M 
d?P 



C 



M 



d 3 R P ■ u(R) 
(2^)3 



dxw a (R;xP+) \<t> M {x)\ 2 w b (R;(l-x)P+). 



(44) 



The amplitude <Pm(x) encodes the remaining QCD dy- 
namics. We expect it to be peaked around x = 1/2, 
meaning that the two quarks will carry roughly the same 
amount of momentum. But the width of the distribu- 
tion, since it is formulated in terms of momentum frac- 
tions, could be quite broad in momentum space. Thus 
we cannot use the same argument as for the transverse 
coordinates in order to integrate out this degree of free- 
dom. However, for an exponential parton spectrum we 



have 

w a (R;xP+)w b (R;(l-x)P+) 

~ e -^ + /T e -(i-*)p+/T = 6 -p+/t (45) 

Hence the product of parton distributions is independent 
of x and we can perform the integral over x, which just 
gives the trivial normalization of the wave function from 
There will be corrections to that from momen- 
tum components other than P + which are not additive 
because energy is not conserved (see next subsection). 
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Where we want to take into account wave functions, we 
adopt the asymptotic form of the perturbative pion dis- 
tribution amplitude 



4>m{x) 
as a model. 



30x(l - x) 



(46) 



C. Energy conservation 

Energy conservation is not manifest in the recombina- 
tion approach. Since we are dealing with a 2— >1 pro- 
cess, one of the particles in general needs to be off mass 
shell [35| . This does not pose a problem in the physi- 
cal environment where recombination takes place. Both 
in the quark phase and in the hadronic phase we ex- 
pect interactions with the surrounding medium to occur. 
Since recombination, as described in this work, has been 
simplified to a counting of quantum numbers and mo- 
menta, without real QCD dynamics, we neglect effects 
of these additional interactions that ensure energy con- 
servation. This is tolerable due to the small time scale 
of hadronization and that changes in the parton distri- 
butions w by interactions during this time are negligible. 
We also neglect effects that final state interactions be- 
tween the hadrons could have on the spectra - these are 
expected to contribute mainly in the low momentum do- 
main, Pt < 1 GeV. 

For large momenta P, the energies of the particles are 
dominated by the kinetic energies and not by the masses. 
The light cone wave functions have a transverse width Aj_ 
which is a non-perturbative momentum scale. Therefore 
the momenta of the recombining quarks will be collincar 
up to transverse momenta of order Aj_ • This implies that 
the energy is conserved up to terms of order A 2 L /P + and 
m 2 /P + , where m stands for the masses of the participat- 
ing particles. 

Energy conservation is a problem at low transverse mo- 
mentum and can only be overcome if one takes further in- 
teractions between the partons into account. This would 
require a much more sophisticated formulation that in- 
cludes nonperturbative initial- and final-state effects. 



Restricting ourselves to large transverse momentum fi- 
nally permits to use light cone fractions x for the spatial 
momentum P in the CM frame instead of the momentum 
P + in the HLC frame. The relation p + — xP + between 
the quark momentum p and the hadron momentum P in 
light cone coordinates translates to 



xP + 0\x 



M 2 



(47) 



This is a leading order expansion in M 2 /P 2 where M 
is the hadron momentum. At mid rapidity this further 
translates to the simple formula px = xPt for the trans- 
verse momenta. 

D. Low transverse momentum and hadron 
thermodynamics 

It is well known that total hadron yields at RHIC 
can be described very accurately by a purely statistical 
model, using only hadronic properties, such as masses, 
chemical potentials and spin degeneracies [HI- 
The hadron yields are given by the P^-integrated spectra 
and are naturally dominated by particles with less than 
2 GeV/c transverse momentum. The picture of thermal 
hadron production does not necessarily require the exis- 
tence of a parton phase. 

However, one can show that recombination of a ther- 
malized parton phase is consistent with thermal hadron 
production in the limit Pt — ► oo. Although pQCD will 
eventually dominate over both mechanisms at large Pt, 
this nevertheless suggests that the recombination mecha- 
nism connects a thermal parton phase with the observed 
thermal hadron phase. Therefore it is justified to call 
recombination from a thermal parton phase the micro- 
scopic manifestation of statistical hadron production. 



E. Summary of the formalism 

We want to summarize what we have so far. From (I44|) 
the meson spectrum is given by 



J 



(48) 



s o 

where dan measures the volume of the hypersurface £ and R G S. For baryons, the same steps result in the expression 



E WF =Cb J d<7R ~727rp~ ^ VXi w 4^ x i P ) w &(^ x 2P)w c (i?;x 3 P) \(j>B{xiiX 2 ,x 3 ,y 2 
s 



(49) 



r 



a, b and c are the valence partons and 4>b{%i, %2, £3) is coordinates. We use the short notation 
the effective wave function of the baryon in light cone 1 

T>Xi = J dx\ dx2 dx^5(xi + x 2 + £3 — 1) (50) 

b 
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for the integration over three light cone fractions. In- 
spired by the asymptotic form of the light cone distribu- 
tion amplitudes for pions and nucleons we choose 



(51) 
(52) 



<I>m(x) = v30a;(l - x), 
b B (xi,x 2 ,x 3 ) = 12V35xix 2 x 3 . 



These wave functions are broad in momentum space. 
In order to study the effect of the width of the wave func- 
tions on our results it will be interesting to alternatively 
explore the case of narrow wave functions in the spirit of 
Sec. Ill Dl The limiting case are (5-shaped wave functions 



l".wU)| 2 = S(x- - 



(53) 
(54) 



The spectra are then given by 



E 



d 3 P 



c 



M 



da. 



P ■ u(R ) 
' (2tt) 3 



- WQ (i?;^) Wb (i?;^) (55) 



s 



■ a (R;j)wb(R)j)w c (R-,j). (56) 



X w, 



It is an important observation that in the case of expo- 
nential parton distributions, the shape of the wave func- 
tion is almost negligible. We have to be aware that there 
will be corrections to the above equations of order m/|P| 
and Aqcd/|P|, where m is the mass of the hadron or the 
partons, reducing their range of applicability to large P. 

On the other hand the spectrum of hadrons from frag- 
mentation is given by 



E 



,dN h 
d 3 P 



E 



dz 



D a ^ h {z)E a 



dN a 
<PP a 



(57) 



The number of partons can be obtained from the cross 
section via the impact parameter dependent nuclear 
thickness function dN/d 3 P = T A uAu{b)da R /d 3 P. We 
take TauAu(O) = 9A 2 /8nR\ for central collisions. Ra = 
A 1 / 3 1.2 fm is the radius of the nucleus. 



IV. HADRON AND PARTON SPECTRA 

In the following we want to discuss the parton phase 
created in collisions of gold nuclei at RIIIC with y/s = 
200 GeV per nucleon pair. We will then proceed to cal- 
culate the hadron spectra emerging from recombination 
and fragmentation of this parton phase. 



A. Modeling the parton phase 



Assuming longitudinal boost invariance, we fix the hy- 



persurface S by choosing r = 



const, for 



W 1 = (t,x,y,z) 

= (t cosh r], p cos 4>, p sin </>, t sinh if) . (58) 

It is convenient to introduce the space time rapidity rj and 
the radial coordinate p, since the measure for the hyper- 
surface £ then takes the simple form dan = r dr\ p dp d(f> 
and the normal vector is given by 



'(R) = (cosh??, 0, 0,sinhf7). 



(59) 



Using a similar parametrization of the parton momentum 



nt 1 — 



{rtiT coshy, px cos $,pt sin $, tot sinhy) (60) 



with rapidity y and transverse mass tot = y m 2 + P T , we 
obtain p ■ u(R) — tut cosh(ry — y). We remind the reader 
that we will use capitalized variables (P, Pt, M, Mt, etc.) 
for hadrons. 

Now we have to specify the spectrum of partons. As 
already discussed above we assume that the parton spec- 
trum consists of two domains. At large pt, the distribu- 
tion of partons is given by perturbative QCD and follows 
a power law. For the transverse momentum distribution 
at midrapidity for central collisions (b=0) we use the 
parametrization [45ll46| 



dNJ! c 



d 2 prdy 



= K °- 

y =o {1+px/BY 



(61) 



The parameters C, B, and f3 are taken from a leading 
order (LO) pQCD calculation and can be found in |4y| 
for the three light quark flavors and gluons. A constant 
K factor of 1.5 is included to roughly account for higher 
order corrections in a s |47| . The calculation includes 
nuclear shadowing of the parton distributions, but no 
higher twist initial state effects. Higher twist effects, 
like the Cronin effect, will fade like A 1 / 3 Aq CD /Py for 
high transverse momentum |48| . Since we will show that 
fragmentation and pQCD are only dominant for trans- 
verse momenta above 5 GeV/c in the hadron spectrum 
for RHIC, it is safe to omit the Cronin effect. 

Energy loss of partons, resulting in a shift of the trans- 
verse momentum spectrum 0, 0| , is taken into account 
and parameterized as 



Ap T (b,p T ) =e{b)^fpT^-. 

Ha 



(62) 



For central collisions we take (L) — Ra and therefore 
Apt{0,Pt) — e (0) v / pr (we postpone the discussion of 
the impact parameter dependence to Subsection IV. C). 
The choice of (L) = Ra neglects the fact that for the 
strong quenching observed at RHIC energies, jet emis- 
sion becomes a surface effect However, we note 
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Pt (GeV) 

FIG. 1: Spectrum of u and s quarks at hadronization in a 
central Au+Au collision at RHIC. Perturbative partons from 
hard QCD processes with subsequent energy loss (dashed 
lines) and the thermal phase with T — 175 MeV and radial 
flow vt = 0.55c (solid lines) are shown. 



that only the product e(L) as a whole is a parameter. We 
have no ambition here to make a connection to the mi- 
croscopic parameters of jet quenching, therefore we will 
not disentangle e and (L) . This would require a more so- 
phisticated model of the emission geometry. We also do 
not use a radial profile for the emission and the density 
of the medium. This can be found discussed elsewhere 
in the literature P. l2ll 12^ . Nevertheless, our "minimal" 
description of energy loss is quite successful to describe 
the available data on high-Py spectra (Pt > 5 GeV/c) 
for 7T°, and charged hadrons (h + + h~)/2 in cen- 
tral Au+Au collisions. From a fit to these data we find 
60 = e(0) = 0.82 GeV 1/2 for central Au+Au collisions. 
This value corresponds to an average energy loss of 3 
GeV for a 10 GeV parton in a central Au+Au collision. 
The perturbative spectrum for up and strange quarks 
at midrapidity is shown in Fig. ^ For fragmentation of 
pions, kaons, protons and antiprotons we use LO KKP 
fragmentation functions with the scale set to the hadron 
transverse momentum Pp |19| . A fragmentation is calcu- 
lated with the LO fragmentation functions of de Florian, 
Stratmann and Vogelsang [Hoj . 

Besides the perturbative tail of the parton spectrum 
that will turn into hadrons via fragmentation, we assume 
the existence of a spectrum of thermalized partons that 
are recombining at hadronization and dominate at low 
and intermediate values of pr- In this phase we assume 
the effective degrees of freedom to be constituent quarks 
without dynamical gluons. We take the spectrum to be 
exponential with a given temperature T 

w a {R-p) = jae- p - v{R)/T e-i 2 / 2A2 f(p, 0). (63) 



7 a is a fugacity factor for each parton species a. We also 
include longitudinal and radial flow through the velocity 
vector 

v ll (R) — (cosh rjL cosh tjt, sinhf^ cos <f>, 

sinhryx sin <j>, sinh T)l cosh rjr)- (64) 

t)l(R) and T)t(R) are the rapidities of the longitudinal 
and radial flow which still could depend on the space 
time point R G S. For the longitudinal expansion we 
choose a Bjorken scenario where the longitudinal rapidity 
is simply fixed by the space time rapidity 

Vl{R)=V- (65) 

The transverse flow is given by a velocity vt(R) with 
vt = tanh rjr- For practical purposes we will not work 
with a radial profile but assume vt to be independent 
of p and <f>. However, for collisions with finite impact 
parameter b we will later allow a dependence of vt on 
the azimuthal angle cj> in order to describe the measured 
elliptic asymmetry in the spectra. The space-time struc- 
ture of the parton source in l|63[) is given by a transverse 
distribution /(p, (j>) and a wide Gaussian rapidity distri- 
bution with a width A. 

We assume that hadronization occurs at r = 5 fm at a 
temperature T = 175 MeV in the parton phase. This is 
consistent with predictions of the phase transition tem- 
perature at vanishing baryon chemical potential from lat- 
tice QCD |4l| . For the spread of the parton distribution 
in longitudinal direction we choose A = 2. The con- 
stituent quark masses are taken to be 260 MeV for u and 
d quarks and 460 MeV for s quarks. 

The two component model of the parton spectrum with 
an exponential bulk and a power law tail is also predicted 
by parton cascades like VNI/BMS 51], although the 
interactions in that case are purely perturbative. This 
implies that an exponential shape of the spectrum does 
not necessarily mean that the parton system is in thermal 
equilibrium. 

In the region where contributions from recombination 
and fragmentation are of the same size we expect other 
mechanisms to play a role, which interpolate between the 
two pictures. This could include partial recombination 
and higher twist fragmentation. In the absence of a con- 
sistent description of these mechanisms we simply add 
both contributions to the hadron spectrum - recombi- 
nation from the exponential part and the fragmentation 
from the pQCD part - for P T > 2 GeV/c. 

B. Degeneracy factors 

It is not a priori clear from QCD what the degeneracy 
factors Ch for each hadron h are. In principle every quark 
has 3 color and 2 spin degrees or freedom. On could 
argue that 3 quarks of any color and spin can form a 
proton and that quantum numbers can be "corrected" 
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at no cost by the emission of soft gluons. That would 
lead to degeneracy factors C p — 2 x (3 x 2) 3 and C n o = 
(3 x 2) 2 . On the other hand there are no dynamical 
gluons in our picture and it would be consistent to require 
recombining partons to have the right quantum numbers 
at the beginning. 

Surprisingly this is supported from work on recombina- 
tion in pQCD, where the contributions from color octets 
and spin-flip states to the recombination of D mesons 
were found to be small |2(j. Using this assumption, the 
degeneracies are only determined by the degrees of free- 
dom of the hadron, e.g. C p — 2, — 1 etc. These are 
exactly the degeneracies used in the statistical thermal 
model. We will not take into account feeddown from de- 
cays of resonances, except for the A, where the E° is too 
close in mass to be suppressed. Hence we use Ca = 4. 

This is different from [T(il | where we counted A res- 
onances to give nucleons with weight 1. However this 
overestimates the correction from A decays. Neverthe- 
less the degeneracy of 5/3 given in 10] for the proton 
was of the right size due to a mistake in the normaliza- 
tion of the baryon states, that gave an additional factor 
of 1/3!. Therefore the numerical results given in |lol | are 
still valid. 

We should add that due to the small but probably non- 
vanishing color octet and spin flip contributions and due 
to feeddown corrections we expect all degeneracy factors 
to have an error of at least 20%. 



where m-r is the transverse mass of quark a. The factor 
g = 6 in (j66|l is counting the color and spin degenera- 
cies. In the case of central collisions (impact parameter 
b = 0) we can assume that vt and / are independent of 
<j). For simplicity we furthermore assume that / is also 
independent of the radial coordinate p and that the ra- 
dial distribution of the partons is homogeneous up to a 
radius pq: f(p) — Q(po — p). We can then easily per- 
form the 4> and p integrals in i|66[l . If we consider the 
region around y = 0, we can also neglect the Gaussian 
profile function in rj, since r] 2 /A 2 <§C (ttit/T) cosh.77, and 
integrate over this variable analytically. 



The transverse momentum spectrum of partons in the 
thermal phase is then given by 



N. 



in 



(Pprdy 



x I Q 



Pt sinh r/T 




rriT cosh r\T 


T 




T 



(68) 



Central collisions 



For the momentum spectrum of quarks 
dlC f p-u(R) 

E ^fp- =9ia ! <kjR i^r Wa{ R;p) 

s 

we rewrite the exponent in l|63|l as 



(66) 



p ■ v(R) — rax cosh(?7 — y) cosh tjt 

— pt cos(0 — $) sinh tjt (67) 



At = Pq^ is the transverse area of the parton system 
and Iq and K% are modified Bessel functions. Fig. ^ also 
shows the thermal spectrum of up and strange quarks. 
The radial flow velocity vt = 0.55c, fugacities ^ u — Jd — 
1, lu = Id = 0.9, 7 S = 75 = 0.8 and the radius po = 
9 fm were determined by fits of our calculation to the 
measured hadron spectra in central collisions, see below. 



The transverse momentum spectrum of hadrons 
formed by recombination from the thermal parton spec- 
trum can be derived from Eqs. (|48I49|I to be 



N M 



d 2 P T dy 



C M M T ^-^2 lalh I 



Pt sinh t/t 



T 



dx \(/>m(x)\ k M (x,Pr) , 



(69) 



d?P T dy 



= C'bMt^^t 27 Q 7h7 c Iq 

y=o {2ttY 



Pt sinh t\t 



T 



Vxi \(/>b(xi,X2,x 3 )\ k B {xi,P T ) 



(70) 
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for mesons and baryons respectively. We introduced the short notations 
cosh TjT 



k M {x,P T ) = K 1 
k B (xi,P T ) = Ki 



T 

cosh r/x 



x 2 P 2 



(1 - x) 2 P 2 



T 



™2 p2 
x l- r T 



~2 p2 
• L 2 r T 



■2 p2 



Here m a , rrib and m c are the masses of the valence quarks 
and Mt is the transverse mass of the hadron. Fig. [21 
shows the result of our calculation for the spectra of 
pions, protons, antiprotons, kaons, Lambdas, Xis and 
Omegas. See Sec^for discussion. 



D. Peripheral Collisions 

In order to describe peripheral collisions, we have to 
scale the perturbative part of the parton spectrum given 
in I jfiljl by the ratio of thickness functions, or equivalently 
by the number of binary nucleon collisions 

dNr«(b) = ^4S^r rt = w^dNr'- (73) 

-tAuAu(U) -ZVcollfU) 

Our values for the number of collisions N co n as a function 
of impact parameter are listed in Tab. [fl and are close to 
the values used by the PHENIX collaboration [52|. 



(71) 
(72) 

I 

two nuclei 

AT{b) = H^ AT{0) = mw{b) ^ (75) 

A A 

In principle, the radial flow velocity vt is expected to 
vary with impact parameter. However, it turns out that 
the slope of the measured hadron spectra above 2 GeV/c 
is consistent with a constant flow velocity vt = 0.55c 
up to very large impact parameters, so that we fix this 
value for all b. It may be questionable whether partons 
are produced in equilibrium in peripheral collisions. We 
therefore introduce an additional impact parameter de- 
pendent fugacity 7(6) common to all quark flavors. How- 
ever, the measured hadron spectra favor 7(6) = 1 up to 
high values of b. Corrections are only necessary for very 
peripheral collisions. We take 7(6) = 1 for b < 10.5 fm, 
7(11 fm) = 0.7, 7(12 fm) = 0.4 and 7(13 fm) = 0.4. 

E. Elliptic Flow 



b [fm] 





5.5 


7.5 


9 


10 


11 


12 


13 


13.9 


A^coll 


1146 


594 


350 


199 


120 


61.6 


26.0 


10.0 


5.3 



TABLE I: Average number of binary nucleon collisions N co n 
for some values of the impact parameter b in collisions of gold 
nuclei. 

The length and width of the overlap zone of two nu- 
clei with radius Ra, colliding at impact parameter b, are 
1(b) = y/R 2 A - (6/2) 2 and 10(6) = R A ~ b/2. We scale 
the average length entering the energy loss in l|62(l as 
(L) = (1(b) + w(b)) /2. On the other hand the density 
of the hot medium is decreasing with increasing impact 
parameter. We choose the simple ansatz 



e(b) = £ 



1 _ e -(2R A -b)/R A 

1 -e- 2 



(74) 



for the b dependence of the energy loss parameter, which 
describes the data surprisingly well. We refer to 3 for 
more detailed studies of the jet quenching effect. 

For the thermal phase of the parton spectrum, we keep 
the temperature T and the hadronization time r indepen- 
dent of the impact parameter b, but adjust the size of the 
volume according to the profile function f(p). We scale 
the transverse area At of the parton phase at hadroniza- 
tion with the transverse area of the overlap zone of the 



For hadrons from fragmentation we assume that the 
azimuthal anisotropy in peripheral collisions is induced 
by the azimuthal dependence of the energy loss 0, 0, 
I53I IH^j . To determine the coefficient t>2 we generalize 

(E2J) to 

A PT (b,p T ,$) =e V^^-(l-acos2$), (76) 
Ha 

where e(b) is given by Ij74|l . This leads to a spectrum 
d 2 N/dP T with non trivial dependence on $. From that 
we obtain the elliptic flow by applying the definition 



V2(Pt 



f d$ cos 2$ d 2 N/dPl 

< c ° s2 *>- ■ (77) 



The parameter a is given by the collision geometry. It 
has the value 



a 



w(b) - 1(b) 
w(b)+l(b) 



(78) 



for given impact parameter b. 

In the thermal phase, the hydrodynamic expansion 
from an originally anisotropic overlap zone of both nu- 
clei (for b ^ 0) induces an elliptic anisotropy of the par- 
ton spectrum. This is leading to a dependence of the 
transverse flow on the azimuthal angle <fi. To model the 
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FIG. 2: Hadron spectra at midrapidity as a function of transverse momentum Pt for central Au+Au collisions at y/S = 200 
GeV. We show fragmentation (dotted), recombination (dashed) and the sum of both contributions (solid line) versus data from 
PHENIX and STAR (for $, H + + H~ and + fl recombination only). Q + Cl data is minimum bias, therefore the result of a 
calculation with impact parameter b — 10 fm is also shown. All data is preliminary except n°, p and p data from PHENIX. All 
error bars give statistical errors only, except for 7r° data from PHENIX which give the total error. See Sec.EJfor more details. 
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anisotropy in the partem phase we take the azimuthal 
dependence of the transverse rapidity to be 



J*rfo) = »#(l-/(pr)cos20) 



(79) 



Here rfp is the rapidity given by the flow velocity vt = 
0.55c, so that tanh^ = 0.55. The amplitude of the 
anisotropy f(pr) is given by the geometrical anisotropy 



a at low transverse momenta, but faster partons will ex- 
perience the anisotropy in the expansion less than slower 
ones. Therefore we choose an ansatz 



f(PT) = 



{pr/po) 2 ' 

From Eqs. I|66I67|) one obtains |5^, |56 



(80) 



v a 2 {pi 



(cos(2$)) 



/ d<j) cos 2<j) I 2 [p T sinh tj t {(f) /T] K l [m T cosh rj T (0) /T] 
J d(/)Io [p T sinh77 T ((/))/T] K x [m T cosh r] T (<j))/T} 



(81) 



Our assumptions imply v 2 = v 2 = v 2 = v 2 in the ther- 
mal phase, but v 2 = v 2 differ slightly at small pt because 
of the bigger strange quark mass. We determine the pa- 
rameter pq in the parametrization of the parton phase 
from a comparison to the PHENIX measurement of the 
elliptic flow of pions [57]] • We obtain p n = 1.1 GcV/c. 

After having fixed the coefficient v 2 (pr) for each par- 
ton species a we write the azimuthally anisotropic phase 



space distribution for thermal partons at midrapidity as 



wf(R;p) = w a (R;p) [1 + 2v a 2 {p T ) cos2$] 



(82) 



Here w a (R',p) is the phase space distribution without 
anisotropy from Eq. (|63[) . Substituting this into the basic 
recombination formulae 1)481 149fl we obtain 



<(P T ) 



J dx\Mx)\ 2 [v$ {xP T )+v b 2 {{l-x)P T )\ k M (x,P T ) 
J dx \<t> M (x)? [1 + 2vg (xP T ) v\ ((1 - x)Pt)] k M (x,P T ) 



(83) 



v B {P T ) 



J Vx t [tg {xiPt) + v\ (x 2 P T ) + v c 2 (x 3 P T ) + 3z; 2 a ( Xl P T ) v\ (x 2 P T ) v c 2 (x 3 P T )] k B {x u P T ) 

J V Xi |<M^)I 2 [1 + 2 (v% ( Xl P T ) v\ (x 2 P T ) + v% ( Xl P T ) v c 2 (x 3 P T ) + v b {x 2 P T ) v c 2 (x 3 P T ))] k B ( Xi ,P T ) 

I 



(84) 



for the anisotropics in the meson and baryon spectra re- duces to the relations already given before in the litera- 
spectively. Using the (^-function approximation this re- ture [l4l fl5j 



M (p s V^P T )+v\(lP T ) 



i + 2v^\p T )v b 2 c 1 p T y 
V2 T) 1 + 2v«ap T )v b 2 ap T ) + 2v b 2 ap T y 2 ap T ) + 2v c 2 ap T )v%ap T ) ' 



(85) 
(86) 



r 



If we assume one universal partonic v 2 for the recom- 



bining quarks the above expressions simplify to 
2« 2 (iP T ) 



V2,m( P t) 



v 2 ,b{Pt) 



~l + 2v 2 { 2 -P T ) 2 ' 
Zv2(\Pt) + 3v 2 (±P T ) 
l + 6v 2 aP T ) 2 



(87) 



1G 



Since the maximal values of will be of the order of 0.1, 
we can neglect the quadratic and cubic terms and arrive 
at the following simple scaling law, which connects the 
elliptic flow of hadrons v\ to those of the partons v^- 



v%(P T )=nv 2 (~ p T^ 



(89) 



with n being the number of valence quarks and anti- 
quarks contained in hadron h. This scaling law was in- 
deed already found to hold in STAR data on the elliptic 
flow of A and K® down to transverse momenta of about 
500 MeV/c ^(|. This is a very strong support for the 
recombination picture. Apparently a part of the uncer- 
tainty in the recombination mechanism at low Pt, in- 
troduced by the violation of energy conservation, cancels 
after taking the ratios in Eqs. ()83I84[I . The recombina- 
tion formalism seems to give valid results for i>2 down to 
transverse momenta of several hundred MeV/c. 

We combine the contributions to the anisotropic flow 
from recombination and fragmentation by using the rel- 
ative weight r(Pr) for the recombination process 

v 2 (Pt) = r{P T )v2M{P T ) + (1 - r(Pr)) v 2 ,f(P T ). (90) 

r(Pr) is defined as the ratio of the recombination contri- 
bution to the spectrum and the total yield. 



r{P T ) 



dN K /d 2 P T 



(dN R /d 2 P T + dN F /d 2 P T ) ' 



(91) 



F. The statistical thermal model 

In this subsection we give a brief account of the sta- 
tistical model following variant I of 0] ■ For further de- 
tails we refer the reader to the comprehensive literature 

El El El m. 

The hadron spectrum at is supposed to emerge from a 
hypersurface II and has the form 



E 



dNh 
d 3 P 



= / da, 



P ■ v(R) 
' (2tt) 3 



G h (R;P). 



(92) 



We use the same parametrization for the four velocity 
v(R) as in i|64[l . The hypersurface II is determined by 
the condition = tsm = const. The hadronic phase- 
space distribution functions are given by 



G h (R;P) 



ChfsM{r) 



(P-v— urBu— UsSh— /ilifc )/Tsm -J- 1 



, (93) 



for bosons and fermions respectively, r = tsm sinh tyr 
is the radial coordinate and /sm{t) — 0(ro — r) is a 
radial profile function providing a cylindrical shape. Ch 
is the degeneracy factor and Bh, Sh and Ih are baryon 
number, strangeness and third component of the isospin 
for hadron species h. 

Equation i|92|l can be evaluated analogous to (|66[l . We 
note that in the limit Pt — > oo Eqs. I|69I70|) are equivalent 



to i|92|) if the same hypersurface and the same tempera- 
ture and chemical potentials are used. This is an indi- 
cation that recombination from a thermal parton phase 
is the underlying microscopic picture of hadron produc- 
tion in a statistical model. While we will not elaborate 
on this in more detail, we will quote some results of the 
statistical model for hadron ratios and compare with our 
calculation. 

The geometric parameters are fixed to be tsm — 7.66 
fm and ro = 6.69 fm for for most central collisions at 
RHIC in Ref. E3- Particle ratios at mid rapidity in a 
boost-invariant model are not influenced by the expan- 
sion of the system E2 , thus we can use the parameters 
which are determined by particle ratios from the entire 
phase space. We follow E3 and set T SM = 177 MeV, 
fi B = 29 MeV, fi s = 10 MeV and /// = -0.5 MeV. 



G. Note on the parameters in our model. 

We want to give a brief summary of all the parameters 
for the parton phase. Essentially we have three degrees of 
freedom for central collisions. Theses are the energy loss 
given by eo (L) , the slope of the exponential part given by 
temperature T and radial flow velocity Vt and the nor- 
malization of the recombination spectrum by the volume 
tAt- In addition there are the parton fugacities. After 
fixing (L) , T and r to physical or at least reasonable val- 
ues, we retain eo, vt and po as true parameters that were 
determined by fitting to the final data given by PHENIX 
for the inclusive ir° spectrum |52j . This is in contrast to 
our previous study where the parameters of the parton 
spectrum were fixed by the peliminary charged hadron 
spectrum [Tof . 

The light quark fugacity was set to 1 in accordance 
with the measured p/ir° ratio and the fugacities for an- 
tiquarks and strange quarks were obtained from other 
ratios. The ratio of fugacities 7u/7u = 0.9 can be trans- 
lated into a baryon chemical potential /is = 27 MeV. For 
other impact parameters, the simple geometric scaling of 
the volume and the number of collisions with b and a 
reasonable ansatz for e(b) describe the data up to b = 10 
fm. Only for very peripheral collision there is the need 
to introduce the new parameter 7(6). 



V. NUMERICAL RESULTS 

In this section we are going to discuss our numerical 
results on hadron production. 

A. Hadron spectra 

In Fig. |3 we show our results for hadron production 
from fragmentation and recombination for impact pa- 
rameter b = in central Au+Au collision at y/s = 200 
GeV. We compare to available experimental data from 
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FIG. 3: Relative deviation r& of the 8 function approximation 
from calculations using wide wave functions for neutral pions 
and protons. 



the PHENIX and STAR collaborations at RHIC. The 
7T° spectrum has been measured by PHENIX up to 10 
GeV/c. The final data were released very recently [§2 
together with final data for neutral pion production in 
p + p collisions up to 14 GeV/c Error bars show the 
total error in the 7r° yield. All our calculations use the ex- 
pressions (I48I49H with realistic light cone wave functions 
(|51|l . For protons and neutral pions we also performed 
the calculations in the 5-function approximation for the 
wave functions (|55I56|I . Fig. [31 shows the relative devi- 
ation rs = (dN — dNs)/dN of the ir° and p spectrum 
in ^-function approximation from the recombination cal- 
culation using realistic wave functions. The deviation is 
less than 22% for protons and less than 12% for pions at 
small Pt- It becomes considerably smaller at larger Pt 
since the violation of energy conservation is less impor- 
tant there. This explains why calculations using I|55l5t)|) 
are often satisfactory. 

Preliminary data on ir + and it~ production are only 
available up to 2 GeV/c 60], but we can expect that the 
global behavior of charged pions is similar to that of neu- 
tral pions. All data except for ir°, p and p shown in Fig. 
13 are preliminary. Only statistical errors are given for all 
sets besided tt°. In the pion spectra we clearly see the 
two Pt domains of hadron production. Above 4-5 GeV/c 
the spectrum is dominated by fragmentation and follows 
a power law. Below 4 GeV/c the spectrum is exponen- 
tial and dominated by recombination from the thermal 
phase. In our calculation the contribution from fragmen- 
tation is artificially cut off below 2 GeV/c (corresponding 
to a lower cut-off of about 4 GeV/c in the parton spec- 
trum) because perturbative QCD loses its validity at low 
Pt- 

In the crossover region the yield is slightly underesti- 
mated, due to our simplified treatment not allowing for 
recombination involving perturbative partons and mixed 
mechanisms We also note that below 2 GeV/c the 
calculated recombination spectrum bends down and un- 



derestimates the data. This effect is caused by neglect- 
ing the large binding energy of pions in our recombina- 
tion formalism, in which pions have an effective mass 
of 2m u ~ 520 MeV compared to the true pion mass 
of 140 MeV. In addition, pions from secondary decays 
of hadronic resonances are an important contribution at 
very low Pt- 

The same effect can be seen in the kaon spectra, where 
we again underestimate the yield for Pt < 2 GeV/c. 
The K® spectrum was measured by STAR up 6 GeV/c 
|6lj . The preliminary K + and K~ spectra up to 2 GeV/c 
are available from PHENIX [6(| , and first results on K~ 
from STAR up to 4 GeV/c were shown recently |6J. The 
K° data above 2 GeV/c can be described very well by 
our calculations. We note that the last four data points 
of the preliminary STAR data on K~ follow a different 
systematics than the rest of the points and also seem to 
deviate from the data. This could indicate a failure 
of proper particle identification in this momentum range. 

Protons and antiprotons have been identified by the 
PHENIX collaboration up to 4.5 GeV/c and the final 
data was published very recently [13 • In contrast to the 
Goldstone bosons considered before, the mass of p and p 
in our constituent quark picture is closer to the physical 
mass of 938 MeV and secondary protons and antipro- 
tons are less abundant. Hence our calculation provides 
a satisfactory description of the spectra even between 1 
and 2 GeV/c. This is true for all hadrons that are not 
Goldstone bosons. The crossover between the recombi- 
nation and the fragmentation process is shifted upward 
for p and p compared with pions, to around 6 GeV/c. 
We alert the reader to the apparent suppression of the 
fragmentation process below 6 GeV/c compared with re- 
combination. This is much more prominent here than in 
the case of pions. In Fig. 0] we show the ratio r(Pr) of 
recombined hadrons to the full calculation, defined in Eq. 
(|91H . The shift of the crossover to higher Pt from pions 
over kaons to protons is obvious. The 50% mark changes 
from 4 to 4.5 to 6 GeV/c. 

Here we need to emphasize an important point re- 
garding the perturbative calculation. The fragmentation 
functions are a non-perturbative input to these calcula- 
tions, and are derived from other experiments. Most data 
about fragmentation functions are from e + e _ annihila- 
tion experiments which do not allow to distinguish be- 
tween quark and antiquark fragmentation. Quarks have 
to be created in pairs, so that only fragmentation func- 
tions like D u+ u^ p +p can be deduced. Additional input 
from semi-exclusive reactions helps to separate the con- 
tributions, but fragmentation functions still require im- 
provement for applications in hadron interactions. The 
KKP parametrization seems to work well for ir° pro- 
duction at RHIC [59| but we anticipate more problems 
for other hadrons, in particular for protons and antipro- 
tons. This suggests a considerably larger theoretical un- 
certainty for the results from fragmentation of all hadrons 
other than pions. For this reason, and due to the lack 
of appropriate fragmentation functions, we do not show 
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FIG. 4: The ratio r(P T ) = R/(R + F) of recombined hadrons 
to the sum of recombination and fragmentation for it (solid), 
Ks (dashed) and p (dotted lines). For protons and pions 
different impact parameters b = 0, 7.5 and 12 fm (from top 
to bottom) are shown. K® is for b = fm only. 

the fragmentation contribution for $, S and fl in our 
calculations. 

Our results for A + A include an equally large contri- 
bution from S° and E°. The preliminary STAR data 
is taken from |6lj . For the multistrange hadron spectra 
S~ + S+, f2 + Cl and <f> - which is supposed to have a 
pure ss valence structure - we present only recombina- 
tion spectra. The strange hadron yields determine our 
value of the strange quark fugacity 7 S = 0.8. Prelimi- 
nary STAR data is available on S~ + S + [6]J and Q + Cl 
(minimum bias) |63|. We also show a calculation for an 
impact parameter of 10 fm which agrees well with the 
minimum bias + data. No spectra for <J> mesons at 
RHIC have been published so far. 

The charged hadron spectrum (h + + h~)/2 is shown 
separately in Fig. [S] pions, protons, antiprotons and 
kaons are taken in into account. This spectrum, includ- 
ing its impact parameter dependence, will be discussed 
in detail in the section on centrality dependence. 

In summary, our calculations using one fixed parame- 
ter set are consistent with all the currently existing data 
from RHIC and allow us to make predictions for future 
measurements. 



B. Hadron ratios 

In Fig. we show hadron ratios. Only statistical errors 
are shown. The systematic errors can be quite large - we 
refer the reader to the cited experimental publications 
for further details. 

One of the main motivations at the onset of our inves- 
tigation was to find an explanation for the surprisingly 
large proton over pion ratios that are of order one above 
1.5 GeV/c. The data on the p/n° and p/ir° ratios stem 
from the PHENIX collaboration 62]. We have already 
shown in a previous publication 10] that recombination 



naturally provides a p/ir° of order one for hadron trans- 
verse momenta up to 4 GeV. In addition we predict that 
a sharp drop beyond 4 GeV/c should be seen when the 
fragmentation process takes over. The value predicted for 
the ratio in the fragmentation domain is about 0.1. At 
small transverse momenta, the statistical model describes 
the data well but continues to rise beyond 4 GeV/c. 

The K + /ir + and K~ /tv~ ratios have been measured 
by BRAHMS @, and the K~/K + ratio is compared 
to data from PHENIX [gjj. For the statistical model 
we probe an additional strangeness fugacity and provide 
curves for values of 1.0, 0.9 and 0.8. For the recombina- 
tion part 7 S = 0.8 was used as discussed above. The in- 
teresting feature of the preliminary p/p data from STAR 
[H and PHENIX [U is that PHENIX implies a flat p/p 
ratio up to 4 GeV/c, while STAR sees a decrease from 2 
GeV/c on. (This may be related to the apparent surplus 
of K~ measured by STAR in the same momentum range, 
indicating a possible failure of charged particle identifi- 
cation beyond 2 GeV/c.) Our calculation predicts a flat 
ratio in this transverse momentum range in agreement 
with the PHENIX data. A linear sum of statistical and 
systematic error is shown for the STAR data. 

Preliminary results for the ratios (A + A)/4_ft'" and 
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FIG. 5: The spectrum of charged hadrons (h + — h~)/2 for 
four different impact parameters 0, 7.5 (divided by 25), 10 
(/200) and 12 fm (/1000) [from top to bottom] in comparison 
with data from STAR and PHENIX in different centrality 
bins. Contributions from fragmentation only (dotted) and the 
sum of recombination and fragmentation (R+F, solid lines) 
are shown. 
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FIG. 6: Hadron ratios p/ir°, p/ir°, p/p, K+/-K+, K-/ir~, K~/K + , {A + A)/4K°, 2(5" + H+)/(A + A) and 2n°/(h + + h') as 
functions of transverse momentum Pt- The calculation for S baryons only takes into account recombination. We show data 
from STAR, PHENIX and BRAHMS and results from the statistical model (dash dotted lines). Where several curves for the 
statistical model are shown these are for different strangeness fugacities 1.0, 0.9 and 0.8 (from top to bottom). 



2(H + 3 + )/(A + A) were also presented from the STAR 
collaboration Our calculations (for S recombina- 

tion only) are in rather good agreement. Recombination 
predicts a large peak in the A/kaon ratio and a sharp 
decrease beyond 4 GeV/c similar to the p/tt° ratio. First 
indications of such a sharp transition can be seen in the 
STAR data. This observation supports the recombina- 
tion picture including a transition to the fragmentation 
regime beyond 4 GeV/c. 

For completeness we also show the ratio of charged 
hadrons to neutral pions reported by PHENIX [63] in 
comparison with our results. Our calculation agrees very 
well with the PHENIX data in the recombination region 



but slightly underestimates the ratio in the fragmentation 
region. 



In summary, our calculations are in good agreement 
with the available RHIC data. The predicted decrease 
in the proton/pion ratio, if confirmed, and first observa- 
tions of a similar drop in the A/kaon ratio will be strong 
arguments in favor of the recombination+fragmentation 
picture. 
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FIG. 7: The spectrum of neutral pions for impact parame- 
ters 0, 5.5 (divided by 5), 7.5 (/25), 9 (/100), 10 (/200), 11 
(/500), 12 (/1000), 13 (/2000) and 13.9 fm (/5000) [from top 
to bottom] in comparison with data from PHENIX in different 
centrality bins. Contributions from fragmentation only (dot- 
ted) and the sum of recombination and fragmentation (R+F, 
solid lines) are shown. 13.9 fm calculation is fragmentation 
(F) only. 



C. Centrality dependence 

Figs. [7| and [S] show the centrality dependence of the tt° 
and Kg spectra. Final results in various centrality bins 
for 7T° have been published by PHENIXp^. The K° 
data are preliminary results from STAR |6l| . The impact 
parameter dependence of the charged hadron spectrum 
(h + + h~)/2 was already shown in Fig. [S] with data from 
STAR and PHENIX [H. 

The impact parameter dependence of the parameters 
in our calculation was fixed by a fit to the 7r° data but it 
is consistent with the kaon and charged hadron data. We 
notice that with increasing impact parameter the hadrons 
from fragmentation come ever closer to the data points 
below the crossover point. Thus fragmentation becomes 
more and more important in peripheral collisions in ac- 
cordance with our expectations. For the most periph- 
eral bin in 7r°, with an impact parameter of 13.9 fm we 
refrained from extracting a recombination contribution. 
In principle the data can be explained by fragmentation 
alone down to 2 GeV/c. We give a recombination contri- 
bution for b — 13 fm although its contribution might be 
disputable already in this centrality bin. 
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FIG. 8: The spectrum of neutral kaons for impact param- 
eters 0, 10 (/25) and 12 fm (/200) [from top to bottom] in 
comparison with data from STAR in different centrality bins. 
Contributions from fragmentation only (dotted) and the sum 
of recombination and fragmentation (R+F, solid lines) are 
shown. 




P T (GeV) 

FIG. 9: The ratio p/n° for three different impact parameters 
0, 7.5 and 13 fm (solid lines, top to bottom) compared to the 
statistical model (dashed dotted line) and PHENIX data. 



Fig. 01 shows the ratio of r(Pr) for three different im- 
pact parameters (0, 7.5 and 12 fm) for protons and pions. 
The systematics confirms that the crossover point shifts 
to smaller Pj- for increasing b and that the weight of frag- 
mentation becomes more and more important. In Fig. El 
we display the impact parameter dependence of the p/ir° 
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ratio. As expected the proton/pion ratio is decreasing 
with increasing b in the recombination region and is un- 
altered where fragmentation is dominating. 

It is an interesting question at which impact parame- 
ter the recombination mechanism becomes negligible for 
transverse momenta above 2 GeV/c. This will happen 
at smaller impact parameter for pions than for protons. 
Related to this issue is the question whether recombina- 
tion contributes to hadron production at central rapidity 
in p + A or d + A reactions, where the produced mat- 
ter is less dense. We still expect recombination to be 
an important mechanism at low Pp. However, there will 
probably be no chemical and thermal equilibration in the 
parton spectrum, and much less flow. The b dependence 
in Au+Au collisions is not a good basis for extrapola- 
tion since flow and equilibration seem to diminish only 
in very peripheral collisions. That makes it difficult to 
give quantitative predictions without data. However we 
would expect that no recombination effects in pion pro- 
duction are visible above 2 GeV/c in d+Au collisions. 



D. Nuclear modification factors 

The nuclear modification factor is defined as the ratio 
of the hadron yield in Au+Au collisions to the one in 
p + p scaled with the number of collisions 



Raa 



d 2 N Au+Au (b)/dP 2 
N coU (b) d 2 N p+p /dP 2 



(94) 



Similarly one can consider scaled ratios of different cen- 
trality bins like central to peripheral 



Rcp = 



N coll {b) d 2 N Au+Au (0)/dP 2 
AUi(O) d 2 N Au+Au (b)/dP 2 ' 



(95) 



In Fig. EH we show the nuclear modification factor for 
neutral pions for two impact parameters, 6 = and b = 
10 fm. We provide ratios taken both with our own p + p 
calculation and with the p + p results from PHENIX [&3] 
and compare to final data from PHENIX . We notice 
that there is an apparent uncertainty in the perturbative 
calculation which makes the two curves using our own 
p + p calculation and the PHENIX p + p results deviate. 
However, both curves are consistent with the data for 
central collisions. The problem is amplified for peripheral 
collisions, where the curve using our p + p calculation 
overestimates Raa below 4 GeV/c. The spread between 
both curves can be interpreted as a typical error to be 
expected in a lowest-order perturbative calculation. The 
nuclear modification factor for pions shows the strong 
jet quenching effect that suppresses the pion yield by 
about a factor of 5 for the highest Pr bins in central 
collision. Recombination predicts a slight increase below 
4 GeV/c which can be observed in the data. However, 
recombination is not able to compensate the loss of pions 
through jet quenching at high Pp. In peripheral collisions 
jet quenching effects are much weaker. 
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FIG. 10: Nuclear modification factor Raa for impact param- 
eters (bottom) and 10 fm (top) . Normalization by p + p 
is via our own calculation (solid) or via PHENIX p + p re- 
sults (dashed lines). Data on Raa are from the PHENIX 
collaboration with point-to-point errors only. 



Fig. II II displays the scaled ratio Rqp for neutral pions 
and protons. The ratio of impact parameters and 12 
fm is used and compared to data from the PHENIX col- 
laboration [6^. The data for protons shows Rcp to be 
between 0.8 and 1 below 4 GeV/c, which is quite sur- 
prising, considering the strong suppression suffered by 
the pions in that momentum domain. As already no- 
ticed above, recombination is more effective for protons 
than for pions. Therefore our calculation for the pro- 
tons yields a similarly value of 0.8 as observed by the ex- 
periment. This implies that protons from recombination 
make up for the loss suffered from jet quenching at inter- 
mediate transverse momenta. Our calculations predicts 
sharp drops in Rqp and Raa for protons and antiprotons 
beyond 4 GeV/c where fragmentation with jet quenching 
start to dominate. 

In Fig. ^] we give Rcp for charged hadrons compared 
to data from STAR [H and PHENIX 66]. The contribu- 
tion from recombination below 4 GeV/c leads to a value 
of about 0.6 which is between the values for protons and 
pions. The observed steep drop to the value attributed 
to jet quenching is well described by our theory. 

Finally in Fig.^|we compiled Rcp for K® and A+A to- 
gether with data from STAR 61] . The different behavior 
of mesons and baryons is again impressively confirmed. 
For the first time experimental data indicate that a steep 
decrease in i?cp for baryons will occur beyond 4 GeV/c. 
The data on A + A suggest a drop to the perturbative 
value even sharper than what our results show. This 
could be due to too less A baryons from fragmentation 
using this particular set of fragmentation functions |50|. 
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FIG. 11: Rcp for neutral pions (bottom) and protons (top) 
given by the ratio of particle yields at impact parameters 
and 12 fm compared to data from PHENIX. 
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FIG. 13: Rcp for K° s (bottom) and A + A (top) given by 
the ratio of particle yields at impact parameters and 12 fm 
compared to data from STAR. 
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FIG. 12: Rcp for charged hadrons given by the ratio of 
particle yields at impact parameters and 12 fm compared 
to data from STAR and PHENIX. 



E. Results on elliptic flow 
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FIG. 14: Elliptic flow V2 in the parton phase as a function 
of transverse momentum pr- The flow in the thermal phase 
(dashed) and a pQCD calculation (dotted line) are shown. 
The solid line interpolates between the two domains. The 
shaded region shows the region where the interpolation takes 
place. 



Fig. 1 141 shows the elliptic flow v 2 {pt) of u quarks before 
hadronization. The contributions from jet quenching and 
from anisotropic flow in the thermal phase are shown 
separately as well. Due to the very different behavior 
in the two domains, v 2 is more sensitive to mechanisms 
which interpolate between the perturbative domain and 
the soft domain. The range of this theoretical uncertainty 
is highlighted by the shaded region in Fig. ^] 

In Fig.E|we provide results for 7r + . The contributions 
from fragmentation and recombination are shown sepa- 
rately. The full calculation interpolates between these 
two curves in the interval between 2 and 4 GeV/c. We 
also compare the result of a calculation using the d- 
function approximation for the wave functions. As ex- 
pected the deviations are small. Our calculations agree 
with PHENIX data on v 2 of tt+ and 7r~ H3- 

Our results on the particle dependence of elliptic flow 



are summarized in Fig. 1161 One can see the different 
behavior of mesons and baryons by comparing protons 
with pions and kaons with As. v 2 for baryons saturates 
at a higher value than for mesons in the recombination 
domain. At higher Pt, when fragmentation takes over, 
the results rapidly approach each other. In our calcu- 
lation, where we do not take into account the binding 
energies, we cannot resolve the splitting between protons 
and mesons coming from the mass difference. Never- 
theless the agreement with data from PHENIX (5?| and 
preliminary data from STAR [7(j is good. In Fig. 1171 
for ch arg ed hadrons is compared with preliminary STAR 
data |73|. This is interesting since for charged hadrons 
the measurements extend up to 7 GeV/c and constrain 
v 2 in the pQCD domain. Note that in this case v 2 was 
extracted by STAR from 4 particle correlations. This is 
supposed to reduce non-flow effects to v 2 [zll in compar- 
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FIG. 15: «2 for positively charged pions. We show recom- 
bination only (dashed), fragmentation only (dotted) and the 
full calculation (solid line). The result of a calculation in the 
5-function approximation (NW) for the wave function is also 
shown (dash-dotted line). Data are taken from the PHENIX 
collaboration |57|. 

ison with the usual reaction plane analysis. 

In Fig.^]we test the scaling law from Eq. (|89|l for pro- 
tons and pions. Protons and pions follow one universal 
curve below 1.5 GeV/c, which is very similar to the flow 
of thermal partons given in Fig. ^] Beyond 1.5 GeV/c 
we predict a transition to the values given by pQCD. The 
scaling law is no longer valid in that domain. 

We would like to emphasize that we expect modifica- 
tions from other hadronization mechanisms in the region 
where we interpolate between the recombination and the 
fragmentation dominated domains. These could be quite 
important in the case of v-i and alter the results in the 
interpolation region, e.g. they could smoothen the transi- 
tion between both domains. Interactions in the hadronic 
phase could alter vi further. 

VI. CONCLUSIONS 

In this work we have presented extensive evidence that 
recombination is the dominant hadronization mechanism 
for central Au+Au collisions at RHIC up to about 4 
GeV/c for pions and 6 GeV/c for protons. We have 
described a covariant framework that permits the cal- 
culation of recombination from a dense thermal parton 
phase using light-cone wave functions for the produced 
hadrons. This formalism is adequate for momenta much 
larger than the non-perturbative scales involved. At 
lower energies, energy and entropy conservation pose a 
serious problem, the solution of which requires a dynam- 
ical, rather than purely kinematic, treatment of the re- 
combination process. We have found that, for practi- 
cal purposes, hadron spectra are well described down to 
transverse momenta of 2 GeV/c for Goldstone bosons 
(•7T, K) and 1 GeV/c for other hadrons. 

Recombination competes with fragmentation from per- 
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FIG. 16: Upper panel: anisotropic flow for p and ir + com- 
pared to PHENIX data [57| . Lower panel: V2 for K + and 
A + A compared to preliminary STAR data (A + A, K°) [fH 
and PHENIX data (K+ , K~) 0. 

turbatively scattered partons. The large energy loss of 
these partons leads to sizable quenching factors, which 
reduce the fragmentation contribution and cause it to be 
buried under soft physics at scales which one would not 
generally attribute to soft physics. However, 4 GeV/c in 
the pion and 6 GeV/c in the proton spectrum correspond 
to a transverse momentum of only 2 GeV/c on average 
for the coalescing partons. 

The interplay of recombination and fragmentation 
leads to interesting effects in particle ratios and nuclear 
modification factors. The proton/pion ratio is naturally 
around one in the recombination regime. For protons 
the unquenching effect by recombination below 4 GeV/c 
is so strong that essentially no nuclear suppression can 
be observed at all in this momentum range. For the pro- 
ton/pion ratio and suppression factors Raa and i?cp we 
expect a sharp drop beyond 4 GeV/c indicating the be- 
ginning of the perturbative regime. 

For the azimuthal asymmetry V2, our calculations de- 
scribe the data well. The different behavior of baryons 
and mesons is above 1 GeV/c can be explained. The 
scaling law (|89|) , derived from the recombination formal- 
ism, is consistent with data up to 1.5 GeV/c. We predict 
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FIG. 17: «2 for charged hadrons. Again we show the con- 
tributions from different mechanisms as in Fig.EHl Data are 
preliminary and taken from the STAR collaboration. 
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FIG. 18: The anisotropy vijn for pions (bottom) and protons 
(top) as a function of transverse momentum pr/n using the 
scaling law (1891 with n — 2 for pions and n = 3 for protons. 
Data points are pions and protons from PHENIX using the 
same scaling law. 



a violation of this scaling law at higher values, coming 
from pcrturbative QCD. 

In this publication we have only considered sin- 
gle hadron production and neglected correlations in 
the hadron emission pattern. The yield of secondary 
hadrons, when triggering on a leading hadron, is a 
promising quantity to provide more information about 
the underlying hadronization mechanism. 

With fragmentation and energy loss alone, no con- 
sistent explanation involving all hadron species can be 
given. In contrast we are able to describe most avail- 
able RHIC data on spectra, ratios, nuclear suppression 
and elliptic flow of hadrons, including their impact pa- 
rameter dependence, for transverse momenta above 1-2 
GeV/c - for vi even down to very low Pt - consistently 
with a very small number of globally adjusted parame- 
ters. As input for the recombination process we use a 
dense phase of partons with temperature T — 175 MeV 
and radial flow velocity vt — 0.55c at hadronization time 
5 fm. All RHIC data shown in this work are consistent 
with the existence of such a phase. 
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